Synchronization and oscillator death in oscillatory media with stirring 
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The effect of stirring in an inhomogeneous oscillatory medium is investigated. We show that the 
stirring rate can control the macroscopic behavior of the system producing collective oscillations 
(synchronization) or complete quenching of the oscillations (oscillator death). We interpret the 
homogenization rate due to mixing as a measure of global coupling and compare the phase diagrams 
of stirred oscillatory media and of populations of globally coupled oscillators. 
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The problem of synchronization of a large population 
of non-linear oscillators has received a great deal of at- 
tention lately due to its applications in a variety of 
physical 0, chemical Q and biological systems 0,O|, so- 
cial phenomena etc. Typically two types of couplings 
are used to describe the interactions in such systems: i.) 
global coupling, where each of the oscillators is coupled 
to all the others, or ii.) local coupling, where only the 
nearest neighbors are interacting. Global coupling is rel- 
evant for oscillators communicating via visual or acoustic 
signals, like flashing fireflies, chirping crickets, clapping 
audiences, and can also be implemented by electric cou- 
pling (see 3] for a recent experiment on synchronization 
of globally coupled electrochemical oscillators). In these 
cases the time necessary for information to spread over 
the whole system is much shorter than the period of the 
oscillations. In general, global coupling leads to synchro- 
nization when the coupling is sufficiently strong and the 
distribution of the natural frequencies is not too broad. 

In a continuous oscillatory medium (e.g. reaction- 
diffusion system) the oscillations at different points of the 
medium interact through molecular diffusion. Since the 
time-scale of diffusive transport on macroscopic length- 
scales (L) is typically much longer than the characteris- 
tic timescale of the oscillations, L 2 /D 3> T osc , diffusion 
is unable to produce synchronized oscillations over the 
whole domain and the only coherent behavior appears in 
form of propagating waves [3, |8( • 

In certain situations the oscillators are embedded into 
a moving medium, e.g. in a fluid flow. Examples are 
oscillatory chemical or biological systems in stirred reac- 
tors (Belousov-Zhabotinsky reaction |8(, metabolic oscil- 
lations in cell suspensions 0), or in geophysical context: 
oceanic plankton populations |{J and chemical reactions 
in the atmosphere [10( transported by large scale geo- 
physical flows. However, the effect of stirring in oscil- 
latory media has not yet been investigated. It is often 
assumed, that strong stirring leads to spatially uniform 
concentrations, and thus the temporal evolution, simply 
described by a set of ordinary differential equations, be- 



comes independent of the stirring process. But in most 
real systems there are inherent inhomogeneities imposed 
by boundary conditions or non-uniformities of certain 
external parameters. This can be due to spatial vari- 
ations of temperature or illumination in a chemical re- 
actor, or non-uniform distribution of sources in environ- 
mental flows. Therefore perfectly uniform concentrations 
are unattainable and stirring effects cannot be ignored. 

In this Letter we investigate the behavior of a stirred 
oscillatory medium, described by a set of reaction- 
advection-diffusion equations 

d t C i + v(r,t)VC i =F i (C 1 ,...,C N ;r) + DAC i , (1) 

where Cj(r, t), i = 1, 2, N, are the concentrations of N 
interacting species advected by an incompressible fluid 
flow v(r, f). The velocity field is assumed to be time- 
dependent, that ensures efficient mixing by the chaotic 
motion of the fluid elements The functions fj de- 

scribe the nonlinear interactions between the components 
(chemical reactions, evolution of biological populations 
etc.) such that the local dynamics has a stable limit cy- 
cle in each point of the medium 

C i =F i (C 1 ,...,C N ;r)^C i (t + T r ;r)=C i (t;r). (2) 

The explicit dependence of the interaction terms on the 
spatial coordinate accounts for inhomogeneities of the 
medium. We consider the simplest form of inhomogene- 
ity, when the medium is composed by identical oscillators 
except for their frequencies, that is non-uniform in space 



F,(Ci, C N ;r) = (1 + 6f(r))F!(C u .., C N ), (3) 

where /(r) and 6 describe the shape and the amplitude 
of the inhomogeneity. We assume, for simplicity, that 
(/(r)) = and (/ 2 (r)) = 1, where ( ) represents averag- 
ing over the domain. 

In the numerical simulations stirring is modelled by a 
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FIG. 1: Time-dependence of the mean concentration (Ci) for 
different stirring rates (v) and 8 = 0.2. 

sinusoidal shear flow with alternating directions 

v x = {A/T) sin(27ry 4- n ), v y = 
for t e [nT, (n + l/2)T) 
^ = 0, v y = (A/T) sin(27rx + <j> n ) 
for t6[(n+l/2)T,(n+l)T) (4) 

(A = 1.4) advecting chaotically all fluid particles within 
the unit square with doubly periodic boundary condi- 
tions. The phases <p n are chosen at random in each half 
period. This ensures that there are no transport barriers 
and all fluid elements can approach any other fluid ele- 
ment in the domain. We note that chaotic advection is 
a generic feature of simple time-dependent flows, there- 
fore the results are expected to be characteristic to a 
broad range of fluid flows. We define the stirring rate as 
v = 1/T, that is controlled by the period of the flow. 

For the oscillatory dynamics the well known Lengyel- 
Epstein model of the chlorine-iodine-malonic acid reac- 
tion (CDIMA) ji| is considered 

F{(Ci, C 2 ) = 1 - Cx - ACxC 2 /{b + Cl) 
F^(C U C 2 ) = a(C x - CiC 2 /(6 + Cf)). (5) 

The chemical dynamics has a uniform steady state, 
(Cl = 0.2, C| = b + 0.04) which is unstable for the pa- 
rameter values used, b — 0.005, a = 0.375, and the only 
attractor is a limit cycle. The shape of the inhomogeneity 
is chosen to be f(x,y) = y/2sm(2n(x — y)). The system 
(3-5) is investigated for different stirring rates v and de- 
grees of inhomogeneity 5. (D = 10 in all simulations.) 

Let us first discuss some special cases. For non-reactive 
components (F(C\, ..,Cjv) = 0) an initially non-uniform 
concentration is homogenized by mixing. In flows with 
chaotic advection the decay of the spatial fluctuations is 
exponential in time (l2T | 

((C-(C)f)~exp(-2at), (6) 

and for long times the spatial structure is dominated 
by the eigenmode of the advection-diffusion operator 
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FIG. 2: Snapshots of the concentration field Ci(x,y) in the 
synchronized (a.)v = 0.4) and the non-synchronized (b.)v = 
0.167) regimes (8 = 0.2) 

C m i X = DA — v • V with the largest (least negative) 
eigenvalue. 

C(r,i) -> (C(r,0)) +e- at $i, C mix ^ = -a,-**. (7) 

Strictly speaking the eigenmodes have a temporal depen- 
dence following the evolution of the velocity field, but 
they are stationary at least in a statistical sense. 

In a uniform oscillatory medium, 5 = 0, the advection- 
rcaction-diffusion problem has a spatially uniform oscil- 
latory solution. This is stable to spatially non-uniform 
perturbations, and numerical simulations suggest that 
it is globally attracting. The decay of the spatial fluc- 
tuations is controlled by the chaotic mixing, as indi- 
cated by the exponential decay of the variance with the 
same exponent as for the non-reactive case. The only 
difference is that there are oscillations superposed due 
to the oscillatory nature of the chemical dynamics, i.e. 
((C - (C)) 2 ) - h(t) exp(-2at) where h(t) is periodic 
with the period of the oscillations of the mean field. 

Let us now consider a fixed amplitude of the inho- 
mogeneity and vary the stirring rate. When stirring is 
strong the mean concentration has an oscillatory time- 
dependence indicating synchronization (Fig.l), but un- 
like in the case of the homogeneous medium, the spatial 
fluctuations do not disappear completely, since there is 
no spatially uniform oscillatory solution for 8 > 0. For 
very fast stirring the spatial fluctuations are weak and 
the oscillations of the mean concentrations are almost 
the same as for a uniform medium. The amplitude of the 
oscillations of the mean field decreases with the stirring 
rate, while the frequency remains the same. On longer 
time scales, there is also a weak irregular modulation of 
the amplitude, that becomes more pronounced for slow 
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FIG. 3: Projection of the concentration field Ci(x,y) on the 
chemical phase plane Ci — C2 in the synchronized (a.) and 
non-synchronized case, S — 0.2. The limit cycle for a single 
oscillator is also shown for reference. 

stirring. When the stirring rate falls below a certain crit- 
ical value, the synchronized oscillations disappear, and 
the mean concentration is almost constant apart from 
small irregular fluctuations. 

Two pairs of snapshots of the concentration field 
C±(x, y), for the synchronized (a.) and unsynchronized 
(b.) case, are shown in Fig. 2. The complex spatial 
structure, characteristic to chaotic mixing, is combined 
with a coherent time-evolution for supercritical stirring, 
while in the slow stirring case the snapshots correspond- 
ing to different times are statistically equivalent. The 
difference between the two regimes is also clearly visi- 
ble on the projections of the concentration field onto the 
chemical phaseplane C\ — C2 (Fig. 3). When the stirring 
is strong the projected concentration fields appear as a 
small clump moving around the limit cycle of a single os- 
cillator, while in the unsynchronized regime they extend 
over a large domain, that remains almost unchanged in 
time. 

To characterize the degree of synchronization we cal- 
culate the standard deviation in time of the spatially av- 
eraged concentrations 




for different stirring rates (Fig. 4). Large values of R in- 
dicate synchronization. R{v 1 5) can be normalized by di- 
viding it with the same quantity obtained for the uniform 
oscillating medium, R = R(5 = 0) (0 < R/R < I). The 
order parameter R/Rq increases sharply above a critical 
stirring rate. In the unsynchronized regime, below the 
critical stirring rate, the order parameter is small. We 
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FIG. 4: The normalized order parameter as a function of 
stirring rate for different degrees of non-uniformity of the 
medium. 

believe, that in this regime R tends to zero in the D — > 
limit, that is analogous to the limit of infinitely many 
oscillators in the case of global coupling. 

The critical stirring rate depends on 8, faster stirring 
is needed for synchronization when the medium has a 
larger spread of the local frequencies. When 6 is suf- 
ficiently large a new regime appears between the syn- 
chronized and unsynchronized states. At intermediate 
stirring rates the oscillations of the mean field disappear 
completely (R=0) and the concentrations become uni- 
form in space. This 'oscillator death' state corresponds 
to the unstable equilibrium of the homogeneous chemi- 
cal system, C* . Thus, stirring in the presence of inho- 
mogeneity can stabilize the unstable steady state of the 
chemical dynamics and suppress the oscillations. This re- 
sults from the competition between the inhomogeneity of 
the medium generating non-uniform concentrations and 
mixing that tends to reduce the spatial fluctuations. 

Similar regimes have been observed in ensembles of 
globally coupled oscillators 0, Q] . In the stirred media 
due to the relative displacement of different parts of the 
medium the neighborhood of each point is changing in 
time, thus regions that are initially far from each other 
can interact at a later time. This defines a characteristic 
timescale of mixing as the time needed to bring pairs of 
points, initially separated by a distance comparable with 
the size of the domain, sufficiently close so that they can 
interact by diffusion. When this time is much shorter 
than the period of the oscillations, there is an effective 
global interaction in the moving medium. Thus the char- 
acter and strength of the coupling can be controlled by 
the stirring rate. An analogous situation occurs in an en- 
semble of oscillators coupled through a network in which 
each node has a small number of connections that are 
changing in time in a random fashion |15| . 

The dynamics of globally coupled oscillators in the con- 
tinuum limit is described by 

dtd = (1 + 5f(r))F t (C u .., C N ) + a g ((Ci) - Q). (9) 

Similarly to the mixed system, in the absence of os- 
cillations {Fi(C\, .., Cn) = 0), global coupling leads to 
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FIG. 5: Phase diagram of the oscillatory medium with stir- 
ring (a.) and global coupling (b.) - the symbols represent: 
synchronization (■), oscillator death (*) and no synchroniza- 
tion (0). The dashed line is the phase boundary based on 
equation (10). 



an exponential decay of the concentration fluctuations, 
Var(C) ~ exp(— 2a g t). Based on this analogy, we inter- 
pret the homogenization rate in the stirred system, a, 
as a measure of the coupling strength resulting from the 
combined effects of advection and diffusion. 

In Fig. 5 we present phase diagrams, both for the case 
of mixing and global coupling. The homogenization rate, 
a, corresponding to different stirring rates have been ob- 
tained numerically by measuring the decay rate of the 
variance of the concentration field for a non-reactive com- 
ponent. In the fast stirring limit (vL 2 /D ^S> 1) the ho- 
mogenization rate tends to be proportional with the stir- 
ring rate. The two phase diagrams have a qualitatively 
similar structure. Both strong coupling and fast stirring 
leads to synchronization, or oscillator death when the in- 
homogeneity of the medium is strong. Slow stirring is 
analogous to weak coupling as shown by the lack of syn- 
chronization in this regime. 

For globally coupled oscillators an approximation for 
the boundary between the synchronization and oscillator 
death phase has been found recently in 0] 



a c {5)=\ {cu/Xy-l S 2 , 



(10) 



where A and u) are the real and complex parts of the 
eigenvalue of the unstable steady state C* . The above 
result is valid when u> ^> A, in our case A = 0.222 and 
u> = 2.878. We find that in the case of global coupling the 
boundary (|10|) agrees well with the numerical results. For 
the mixed system, however, the same boundary is shifted 
toward lower stirring rates/stronger inhomogeneity. 



This is not surprising since the advection-diffusion op- 
erator has a complex structure and is not simply equiv- 
alent to a linear relaxation to the mean concentration. 
Although the long time decay of the fluctuations is dom- 
inated by the eigenmode with the largest eigenvalue, 
a = ax, by replacing mixing with global coupling cor- 
responding to the most slowly decaying eigenmode we 
underestimate the strength of the coupling due to mix- 
ing. This may be the explanation for the difference be- 
tween the two diagrams in Fig. 5. Another origin of the 
deviation could be that in the system with mixing fluid 
parcels move during a period of the flow and therefore 
the average value of the shape function f(x,y) calcu- 
lated along fluid trajectories has a smaller variance than 
in a motionless medium, resulting in a weaker effective 
inhomogeneity of the medium. 

In summary, we have shown that oscillatory media 
with stirring exhibit qualitatively similar behavior to 
populations of globally coupled oscillators, and the effec- 
tive coupling strength is controlled by the stirring rate. 
Changes in the stirring rate can lead to transitions to syn- 
chronization or oscillator death. This may explain some 
of the stirring effects observed in laboratory experiments 
[lfij | and could also be exploited for controlling the dy- 
namics of oscillatory systems. Similar behavior may also 
arise in systems where the local dynamics is chaotic. 

We thank Silvia De Monte, Peter H. Haynes and 
Francesco d'Ovidio for useful discussions. 
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